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Abstract 

We develop and illustrate methods to compute all single particle potentials that underlie the logistic map, x ^^ sx {1 — x) 
for < s < 4. We show that the switchback potentials can be obtained from the primary potential through functional 
transformations. We are thereby able to produce the various branches of the corresponding analytic potential functions, which 
have an infinite number of branch points for generic s > 2. We illustrate the methods numerically for the cases s = 5/2 and 
s = 10/3. 

I. INTRODUCTION 

In two previous papers [l| it was shown how functions defined on a discrete lattice of time points may be smoothly 
interpolated in t, for a continuum of time points, through the use of solutions to Schroder's nonlinear functional 
equation ^j . If the effect of the first discrete time step is given as the map 

O' 

CN ■ xi^ fi{x,s) , (1) 

^Ij*' for some parameter s, then Schroder's functional equation is 

"^ s*(x,s) = *(/i(a;,s),s) , (2) 

C^ I with ^ to be determined. So, /i (x, s) = ^^^ (s'i' {x, s) , s). A continuous interpolation between the integer lattice 
^ . of time points is then, for any t, 

^' ftix,s) = ^-'{s''^{x,s),s) . (3) 

Oh- . 

This can be a well-behaved, analytic and single- valued function of both x and t provided that ^ ^ (x, s) is a well- 
behaved, analytic, single-valued function of x, even though ^ (x, s) might be, and typically is, multi-valued. In this 

Cd i sense, analyticity in x leads to analyticity in t. 

As discussed in [l[, the interpolation can be envisioned as the trajectory of a particle passing through the initial x, 

x[t)^ft{x,s) , (4) 

^ ' where the particle is moving according to Hamiltonian dynamics under the influence of a potential, V . Up to additive 
^^ , and multiplicative constants, at various times during the evolution of the particle, we have 

^^ ' At t = this becomes V {x) = — (It? s) ( d^(xs)/dx ) ■ Thus, V may inherit multi-valuedness from VP [l[. 

At other times the x dependence of the potential also follows from that of the velocity profile of the interpolation, 
dx (t) /dt, when the latter is expressed as a function of x (t). In general, this will exhibit the branches of the underlying 
analytic potential function. But more importantly for our purposes here, the various branches of the potential can 
also be determined directly from the functional equation V inherits from ^. This functional equation is 



^ 



X 



V{fi{x,s),s)=(-^fi{x,s)] V{x,s) . (6) 

If the map ([T]) possesses a fixed point, we may attempt to solve this functional equation for V by series in x about 
that fixed point. We shall discuss in some detail the circumstances for which this series method is successful in the 
context of the logistic map x >—>■ sx {1 — x) . In general, if the series can be constructed, it will of course have a finite 
radius of convergence. However, the series result can then be continued to other x by making use of the functional 
equation itself (a technique very familiar, e.g., for the F and ( functions) and also by exploiting other special features 
for specific maps (cf. s — >■ 2 — s duality for the logistic map, discussed in Appendix A). These additional techniques 
will allow us to construct convergent series approximations for all branches of the potential in those situations where 
V is multi-valued. The net result is a family of potential sequences that encode for the corresponding continuous 
particle trajectories the various fixed points, bifurcations, limit cycles, and chaotic behavior of the discrete logistic 
map, for all s of interest. 
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II. FUNCTIONAL METHODS AND SERIES SOLUTIONS 

Consider in detail the logistic map [S-Q on the unit interval, x € [0,1] 

X t-^ sx {1 — x) . 



(7) 



For the most part, we restrict our considerations to parameter values s € [0,4]. The maximum of the map is s/4, 
obtained from x = 1/2, so without loss of any essential features, we will also most often restrict x G [0,s/4]. The 
map has fixed points at a; = and at x* = 1 — 1/s. Schroder's equation for this map is 

s* (a;, s) = * (sa; (1 - a;) , s) , (8) 

and from this follows the functional equation for the underlying potential, 

V {sx (l-x), s) = s^ (1 - 2xf V (a;, s) . (9) 

Applying the method of series solution about x ~ 0, with initial conditions that correspond to those used in [Ij for 
the function ^ = a; H — ■' 



+ • • • , namely, V (0, s) = 0, V (0, s) = 0, and V" (0, s) = -2 In s, we find: 



V (x, s) = - (In^ s) U {x, s) , 

/ OO 

U{x,s) ^ x'^ il + y^ an (s) x" 



5 -3s 



n=l 



ai = 



l-s 
The higher coefficients in the expansion are determined recursively by 

n+l 



02 = 



{s -!)'{, + I) 



(10) 
(11) 



1 



hi+2 



(1 - s"+2) 



4a 



n+l 



4a„+ J2 (-ir'«.- 



J + 2 
n + 2 - j. 



for n > 1, 



(12) 



where [• • ■ J is the floor function. In principle, this series solves (El for any s, within its radius of convergence. 
Based on numerical studies, we infer the radius of convergence [8| of the series depends on s as follows: 



i?(s) = 



lim sup 



{\an(s)\^ 



/n 





y 


1 
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= < 




1-i 

s 




s 
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if < s < I 



if I < s < 2 



if 2 < s < 4 . 



(13) 



For < s < 2/3, and also for 2 < s < 4, the |a„| are monotonic for large n and the radius immediately follows 
either from the limsup expression in ([T^ or from the simple ratio test, R{s) = lim |a„_i (s)/a„(s)|. But for 

n— foo 

2/3 < s < 2, there is spiky behavior in |a„_i/a„| for intervals in n, which makes it difficult to use the simple ratio 
test to deternfine R. This is because the |a„ (s)| are not monotonic functions of n for these values of s. Occasionally 
the coefficients become small before changing sign. Fortunately, the lim sup expression for R circumvents this spiky 
behavior to yield the values given in (J13p for all s. 

It is not difficult to work out explicit series results for U {x, s), for generic s, say to O (a;^^). Such polynomials 
in X are sufficient approximations to obtain the graphics to follow, when augmented with functional methods to be 
described. Based on those explicit results, we infer that the series involve numerator polynomials, p„ (s), of order 
1 + (n — 2) (n — 1) /2 in s, as well as "s-factorials" in the following form: 



Pn (s) 



{s-lfU{x,s)^x^is-lf-2is-l)x + J2Yl 

V n=2 ^""^ 



(14) 



Here, deformed integers and factorials are defined by: [k]^ 



-J-, and [n\J = Y[k=i IMs- "^^^ recursion relation for 



the polynomials follows from that for a„ (s). It involves a mix of ordinary and deformed integers: 
1 I _, ,..,.. '^ (-1)""' 



Pn+2 (s) = 



1 



Apn+1 (s) - 4 [n + 1]^ Pn is) + [n + l]J 



^=i+L^J 



b1. 



J + 2 
n + 2-j 



Pj (s) 



(15) 



This recursion relation is seeded by 

Pi(s) = 2(l-s) , p2(s) = 5-3s. (16) 

As written, it looks rather miraculous that the j^ prefactor in ([T5|) is always canceled. Nevertheless, it is. This 
follows from limg^i p„ (s), but we have not yet determined an elegant proof of this fact. 

As originally obtained by Schroder, there are three closed-form solutions known for vjf, for s — —2, 2, and 4. These 
are: 

* (.T, -2) = -^(27T-3 arccos U - l)) . *~' i^^ ^'^^ = \' '=°'' (5 ^ ^ 

* {x, 2) = -1 In (1 - 2x) , VI/-1 (:,, 2) = 1 (l - g-^^) , 

* (a;, 4) = (arcsin Va;) , *""^ (a;,4) = (sin v^) . (17) 

Note that while the VE* are multi-valued the inverse functions are all single-valued. The corresponding closed-form 
expressions for 

U{x,s)= i ,,^}''''\, \ (18) 

are also multi- valued and follow immediately: 

t/(a;,s = -2) = — (l + 2a;)(3-2x) (27r-3arccos(a;- - j j , (19) 

C/ (a;, s = 2) = i (1 - 2xf In^ (1 - 2x) , (20) 

U {x,s ^ 4) = x{l ~ x) arcsin" y/x . (21) 

Were we to start from these expressions for [/, we could recover ^ by solving and integrating (|18p. Another way to 
express the result for s = 4 is similar in form to that for s = — 2, namely, 

1 9 

C/ (a;, s = 4) = -X (1 - x) (tt - arccos (2x - 1)) . (22) 

Indeed, it is well-known that the logistic maps for s = 4 and s — —2 are intimately related through the functional 
conjugacy of the underlying Schroder equations. (See Appendix A.) 

III. OBTAINING THE SWITCHBACK POTENTIALS 

The sequence of switchback potentials, i.e. the various branches of the analytic potential function, can be obtained 
from the functional equation for the potential. (The same procedure also works to give the branch structure of the 
Schroder ^ function. For an example, see Appendix B.) This follows from ([5]) and ([T5)) . namely, 

U {sx (1 - x)) = s^ (1 - 2xf U (x) . (23) 



Next we write y = sx {1 — x), so then x± = ^ I 1 ± y/l — 4j//s j and (1 — 2x±) = 1 — 4y/s. Thus the previous 
relation is just C/ (y) = s (s — Ay) U {x±). Now, we rename y — > a; to obtain 



1 , 1 



x± = -±-^/l-Ax/s, (24) 



U± {x) = s {s - 4:x) U {x±) :^ s {s - 4:x) U (^±^^/l - Ax/s) 



(25) 



One of these potentials (U-) reproduces the original series for the primary potential expanded about a; = 0, as may be 
seen by direct comparison of the series or by numerical evaluation. Alternatively, the other (C/+) gives the potential 
on another sheet of the function's Riemann surface, hence the first switchback potential for s > 2. This is easily 
checked against the closed-form results for the s = 4 case. 

Since C/_ has built into it zeroes at both a; = and a;max = s/4, it is actually a more useful form than just 
the series about a; = 0. When 2 < s < 4, the series has radius of convergence R (s) = s/4. However, from 



using U- instead of the direct series, we have convergence over the whole closed interval, x G [0, s/4], since then 
< ^ — i-\/l — Ax/s < 1/2 < R{s) when 2 < s < 4. Thus we need only evaluate U appearing in U- within the 
region of convergence of its series about zero. That is to say, first construct the series for U, and then from that 
series build C/_. Finally, identify this with C/q, the primary potential in the sequence. 

Similarly, C/_|_ may be identified with the first switchback in the sequence, Ui, but for better convergence properties 
it is useful to build U+ from Uq instead of U. By doing this when 2 < s < 4, wc have convergence over the whole 
closed sub-interval, x G [jss"^ (4 — s) , s/4] , with zeroes of Ui built-in at the end-points of the interval. 

So the primary potential and the first switchback are well-represented by 

f/o {x) = s{s-Ax) u(^^- iv/l-4x/,s^ , (26) 

f/i (x) = s (s - 4x) Uq Q + ]^^l-Ax/^ . (27) 



In the first of these expressions, U is given by the direct series solution of the functional equation 

Now, go through this procedure all over again, beginning with U {sy (1 — y)) = s^ (1 — 2j/) U{y). Let 
z ^ sy{l- y). Then y± = i fl ± y/l ~ Az/s) and U (z) = s {s - 4z) U (y±) = s {s ~ 4z) s {s ~ 4y±) U ((a;±)±) 

where x±± = 2 ( 1 ^ y^l — 4:y±/s ) with any combination of ±s allowed. Therefore U (z) = 



s{s ~ 4z) sfs — 2fl± y/l — 4z/sj j X C/(i±iwl — |fl± y/l — Az/sj I . Again, rename z — ^ a; to obtain 

U±± {x) = s{s- 4x) U± {x±) ^s{s- Ax) s (s - 2 (l ± ^1-Ax/s^'^ ul^±^^l-^[l± ^1 - Ax/s^ J , (28) 

with any combination of ±s allowed. And to improve the convergence properties of this expression, replace U on the 
RHS with Uq. 

This process may be continued indefinitely, through successive application of the basic substitution U (x) — ^ 
s (s — 4a;) U{x±). For example, the next set of potentials in the sequence is U±±±{x) = s {s — Ax)U±±{x±), 
etc. In general, the nth iteration of the procedure gives 

U±± . . .±{x) ^ Hs - ^x)U±± , , ,±{x±) . (29) 

^ ^ ^ ^ ^ ^ 

Finally, at each iteration, we must select appropriate switchback potentials out of the 2" different expressions. In 
particular, we note that many of the U±±,,,± (x) will be complex-valued for the x intervals under consideration, and 
therefore they are not of immediate interest since they do not govern the particle's evolution along the real axis. (The 
continuation of the particle trajectory into the complex plane is outside the scope of this paper.) 

IV. NUMERICAL EXAMPLES 

A. The potential sequence for 2 < s < 3 

These are values of the parameter for which the discrete logistic map converges to a single fixed point: There are 
no bifurcations. Nevertheless, there are two sign choices when the functional equation for the potential is applied 
once, four choices when it is applied twice, and so on. 



V± {x, s)=s{s-Ax)vA\± J] - -, s 1 (30) 



2 V 4 s 



V±± (x, s) = s (s - 4a;) s (s - 2 (l ± ^l-Ax/.^) vA]^± -^--i 



(31) 



As it turns out from numerical studies, for 2 < s < 3 only the positive sign choices are needed to produce the sequence 
of switchback potentials. 



V,{x,s)^s{s~Ax)vA\ + ^j\--^.s\ , v^{x,s)^s{s~Ax)V,r- + ^j\~-^,s\ , (32) 



etc. In general, there is a recursion relation, 



14+1 {x, s) = s{s- Ax) K 2 + 



i-7'^ 



(33) 



The evolving particle moves through this sequence of potentials in succession, with the potential index incrementing 
up by one each time the particle encounters a turning point. 

Consider the specific case s = 5/2. This is representative for 2 < s < 3. For this case, the turning points converge 
onto the nontrivial fixed point x, = 1 — 1/s = 3/5. This is evident in the following graphs. 
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Vq in red, Vi in green, and V2 in blue, for s — 5/2. The first upper turning point is a; = 5/8. Subsequent lower and 
upper turning points are obtained just by iterating the s = 5/2 logistic map, starting with x = 5/8, and are shown 

on the next magnified graph as colored points. 
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Vq in red, Vi in green, V2 in blue, and V3 in orange, for s = 5/2. The first upper turning point is a; = 5/8, the first 

lower turning point is x = j^ ~ 0.58594, the second upper turning point is a; = jfTel ~ 0.60654, etc., as obtained 

by map iteration. The nontrivial fixed point is at x* = 3/5, indicated by the black circle on the x axis. 

As the zero-energy particle moves through this sequence of increasingly shallow, narrowing potentials, its average 
speed decreases, giving the appearance of a dissipative system. Nevertheless, even as the particle motion subsides 
upon convergence into the fixed point at a; = 3/5, energy is rigorously conserved through changes in the potential. 

Insofar as the turning points are branch points for the corresponding analytic potential function, and the various 
switchback potentials are just the values of that analytic function on the various sheets of its Riemann surface, this 
is convincing numerical evidence for that function to have an infinite number of such branch points, for generic s. 



While of course there are many weU-known functions with this property (for example, the inverses of Bessel functions, 
J~^ (x)), this would seem to account for the historical fact that closed- form solutions have not been found for generic 
s, or even for specific s, except in those very special cases where the number of branch points is one or two. Recall the 
exact closed-form cases s — 2 and s — A have analytic potential functions with one and two branch points, respectively, 
as evident in ^^ and (PT|). 

We may check numerically the transit times for each potential in the sequence: Ai„ = 1, for n > 1, within the 
expected uncertainties for truncation of the series (|lip . For example, using Mathematica to compute the initial series 
([TT|) to 200th order and also to evaluate the numerical integrations, we find ^^j^ /^ = 1.000000, /, dx _ 



1.000000, /;5%f^« 

fx(l-a;)l - ^ 

etc 



dx 



1.000000, /^ 

x) 



19 875/32 768 
281241875/2147483648 
J5. 5 /i_ ^| 

128 ' 2 y )\ 



-Vs, 



1.000000, etc. 

19875 5^(i_ 



Vo -•"""""", J75/128 y^VT 

corresponding to the iterations 

_ 1281241875 
19 875/32 768 



X 



1/2 8' 2-^v-^ '^^l2;=5/8 128' 2-^V-^ "'^ 1 2;=75/128 32 768' 2-^">-" ■^'' |2; = 19 875/32 768 2147 483648' 

These transit times are consistent with the continuously evolving, zero-energy particle moving under the influence 
of each of the potentials in succession, Vq ^^ Vi ^ V2 ^ ■ ■ ■ ^ thereby converging onto the fixed point at x = 3/5. The 
numbers confirm that for cases with 2 < s < 3, only positive roots are needed to obtain all the switchback potentials, 
as previously indicated in (|55| . 



B. The potential sequence family for 3 < s < 4 



These are values of the parameter for which the discrete logistic map produces bifurcations and asymptotes to limit 
cycles, rather than unique fixed points, and as s is increased, chaotic behavior erupts. Generally speaking, for these 
values of s the analytic potentials driving the trajectories are more complicated than for s < 3, and their branches 
are not so easy to enumerate as they are encountered by the evolving particle. The potentials here may be grouped 
into a family of sequences^ where potentials with common upper and lower turning points constitute a single member 
of the family. From one family member to the next, for fixed sequence index, the potentials exhibit a progressive 
shallowing (similar to the behavior of the individual potentials for s = 5/2). However, the Vs for a given family 
member successively deepen as the sequence index increases. Consequently, the locally averaged speed of the particle 
moving through the potentials does not necessarily diminish, but often increases, giving the appearance of a driven 
system. Nevertheless, energy is still rigorously conserved. 

The limit cycle situation may be illustrated by the case s = 10/3, while the chaotic situation is aptly illustrated 
by the uppermost value, s — A. The latter case was fully discussed in [1|, and originally led to the notion of the 
switchback potentials. On the other hand, the potential structure for the s — 10/3 case is more complicated than 
for s — 4:, the latter having only one member for its family of potential sequences, and in fact s = 10/3 is much more 
representative of the situation for generic parameter values in the range 3 < s < 4. Therefore, we discuss s = 10/3 in 
some detail. This case evolves towards a two-cycle, consisting of the points ^vTS -I- ^ — 0.830 278 and ^ — ^\/T3 = 
0.469 722 . The first few terms in the potential sequence ((33)) are shown here, for this particular s. 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 




Vo in red, Vi in green, V2 in blue, V3 in orange, and V4 in purple, for s — 10/3. The upper turning point is x = 5/6. 

The two-cycle points are small black circles on the axis. The small black square is the nontrivial fixed point of the 

map, x^ — 7/10, and is revisited after integer time steps, but it is not a fixed point for the continuous evolution. 



Deepening of the Vs is evident for this "mother sequence. 

this s ~ 10/3 case. Rather, they are fixed at | for the upper turning point, and || for the lower turning point, the 



The turning points for this sequence do not converge in 

25 

54 



latter obtained by acting on the former with the s = -y discrete logistic map : | = 0.833 333 i— > If = 0.462 963. 
While these points are not far from the two-cycle points, as evident below, close does not count (except when playing 
petanque) . 
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Vq red, Vi green, V2 blue, V3 orange, & V4 purple, for Vi green, V2 blue, V3 orange, & V4 purple, for s — 10/3. 
s = 10/3. The upper turning point is x = |. The lower turning point is a; = g|. 

Other sequences in the family of potentials are needed to achieve convergence of the continuous particle trajectory, 
for unit time steps, onto the two-cycle of the map. 

So, we define a second sequence in the family of potentials by changing the sign of the root in the arguments of 

the potentials in the first sequence, while noting that Wq (x, s) = s {s — 4x) Vo [^ — \/ j — ^,s\ just reproduces Vq. 

More interesting cases of these "offspring" potentials are given by 



(34) 




These do not reproduce any of the previous potentials. We plot a few for the s — 10/3 case. 
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Wi green, W2 blue, W3 orange, and W4 purple, for s = 10/3. 

For this second sequence the lower turning point can be obtained by acting with the discrete map on the lower turning 
point for the first sequence, i.e. || = 0.462 963 i->- |||| = 0.828 761, while the upper turning point remains the same 
as for the first sequence, i.e. | = 0.833 333. Continuing the enumeration of the potential family members, we define 
a third sequence of potentials, and we plot a few. 



Xn{x,s) = s{s-4:X)Wni- + J- ,S 



for n > 1, 



(35) 
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Xi green, X2 blue, X3 orange, and X^ purple, for s = 10/3. 

For this third sequence the lower turning point remains the same as for the first sequence, i.e. || = 0.462 963, while 
the upper turning point can be obtained by acting with the discrete map on the lower turning point for the second 
sequence, i.e. f§| = 0.828 761 ^ mffff = 0.473 054. 

This procedure may be continued indefinitely. We define fourth and fifth sequences as^ 



y„ (X, s) = S{s- Ax) Xn - 



Zn [X, S) ^ S(s~ Ax)Yn\- + 



1 X ' 

4-7'^ 



1 X ' 

4-7"^ 



(36) 

(37) 



etc., and we plot the first few of these new sequences. 
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Yi green, Y2 blue, and I3 orange, for s = 10/3. 



TA more systematic notation would employ two potential indices, to denote an infinite matrix of potentials. For 

example, we could designate Vn = Vn \ W„ — Vn, Xn — Vri , Yn — Vn , Zn = Vn , etc. The super-index here is 
the family number. 
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Z(x) 
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Z\ green, Z2 blue, and Z3 orange, for s — 10/3. 



The story here is somewhat anafogous to, but more comphcated than, the situation for 2 < s < 3. Instead of 
individual potentials in the sequence {Vn] converging onto the fixed point, as would be the case for 2 < s < 3, here 
each of the individual Vn has proliferated into a "sideways" sequence of potentials, {V„, Wn,Xn, K„, Z„, • • • }, which 
converge to the points in the 2-cycle. To show what we mean by this, we plot {Vi, Wi, Yi} and {Vi, Xi,Zi} near the 
lower and upper points, respectively, in the two-cycle. 
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Wi , and Yi near the lower point in the two-cycle for Vi , Xi , and Zi near the upper point in the two-cycle 
s = 10/3. for s = 10/3. 



This type of behavior is repeated by {Ki, Wn,Yn} and {V^,X„, Z„} for each n. 

To complete the picture for the evolving particle, it is necessary to understand the order in which the potentials 
act, i.e. to determine how the switches are made from one potential to another as the turning points are encountered 
by the particle. Unlike the situation with s < 3, here this can be a bit tedious. From a numerical examination of 
various transit times (and also from the branch structure of the underlying Schroder auxiliary function, as exhibited in 
Appendix B) we infer the order of the potentials for the s — 10/3 case to produce a chemin des energies potentielles: 



Vo^ Vi -^V2^Wi^ W2 
^^^^ ^ f '' ^ 

At=l Ai=l 



¥3^X1^ X2 



V4 



w. 



Yi^Y2^W4,^V5-^ X3-^ Zi^ 



(38) 



At=l 



At=l 



At=l 



etc. Here we have also indicated how the potentials combine into groups with unit total transit time. For example, 
again using Mathematica to compute the initial series (jlip to 20Gth order and computing the transit times numerically. 



we find: J, 



5/6 



dx 



1.000000; /. 
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25/54 y^TT 
dx I /■13575625/28697814 
+ J2J 



5/6 
25/54 



-+I: 



5/6 



dx 



5/6 



V2 ' J 3625/4374 y 

0.164433 



-Wi 



0.825728-1-0.174272= 1.000000; /3g25/4374 



/3W2 



0.661295 + 0.174272 = 1.000000; etc. 



V:, ' J 25/54: ./^X 

It is not difficult to extend (j38p by defining additional sequence family members. Note the turning points must 
match-up for adjacent potentials along the path. Also note the family number plus the sequence number (i.e. 



N — m + n for the potential matrix element Vn of the last footnote) is the same for each potential belonging to a 
group with total At = 1, and this N increments by one as the particle moves from one potential group to the next. 

V. CONCLUSION 

The point of view supported in this paper is that the logistic map, and other discrete time-stepped dynamical 
models, may be regarded as continuously evolving Hamiltonian systems sampled at integer times. For this view to be 
valid, the continuous system must be allowed to undergo a series of switchbacks whereupon the potential affecting the 
dynamics changes when the evolving particle encounters a turning point. From a perspective of configuration space 
covering manifolds, in the case of simple one-dimensional motion, the particle moves from one sheet of a Riemann 
surface to another, to experience a different branch of the underlying analytic potential. 

The methods of this paper may be used directly to determine such branches of the potential for the logistic map, 
for any value of s, as well as other one-dimensional maps. A more extensive study of other examples is underway 
[9|. While peculiar behavior is possible for exceptional maps (say, for special values-!- of s), so far as we are aware, 
any such behavior can always be analyzed using the potential framework presented here. 

In total, for all parameter values governing a particular map, the collection of potential sequence families constitute 
what we may call a "potential fractal" with self-similarities qualitatively evident in the various graphs, as visible 
in the above. Perhaps such potential fractals have a significant role to play in continuum physics. Applications 
might involve any of the usual systems exhibiting chaotic behavior [7|, including accelerator beams jlOj, or perhaps 
cosmological models jll| . 
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VI. APPENDIX A: FUNCTIONAL CONJUGATION 

Functional conjugacy for the logistic map may be expressed in terms of a linear function. Namely, 

g {sx (1 - x)) - (2 - s) g (x) {I - g (x)) , (39) 

where 

g [z) = -^ {l-s + sz) , 5-1 {z)^-{s-l + {2-s)z) . (40) 

Z — s s 

It is often useful to refer to these maps as g (z) = gs {z) and g^^ (z) = g2-s {z)- Pursuing the conjugacy a bit farther, 
let fs (x) ~ sx{l — x), then the conjugacy equation for the map can be written in various ways: 

g° fs ^ .f2-s° g , .fs°g^^ ^ g^^ ° h-s , (41) 

9° Is°g^^ = h-s , g^^ o h-8°9^ Is ■ (42) 

Moreover, Schroder's equation is the functional composition 

s o * = * o /, , (43) 

where s : x — > sx is the simple multiplicative map. Under functional conjugacy by g, the RHS of Schroder's equation 
becomes 

go^o f,og-^ ^go-^og-^o f2-s ^"^9° f2-s (44) 

*gEEgo*og-l (45) 

A. Alternate fixed point expansions 

Consider again the logistic map, x 1— >■ sx (1 — x). On the one hand, with a subscript to distinguish other construc- 
tions to follow, expansions about the trivial fixed point at x = stem from: 

s*o {x, s) = *o (sx (1 - x) , s) , (46) 

$0 (sx, s) = s$o {x, s) (1 - $0 {x, s)) . (47) 

For the inverse Schroder functions, we consistently use the notation $ = ^^^. On the other hand, expansions about 
the non-trivial fixed point at x = 1 — 1/s stem from: 

A*, (z) = *, (Az + (A-2)z2) , (48) 

$, (Az) = A$, (z) + (A - 2) $2 (^) ^ (49) 

where 

1 1 - s 

A = 2-s, s = 2-A, x = l hz, z ^ x^ . (50) 

s s 

The two expansions produce the same functions^ only slightly disguised. Here is the detailed relation between them, 
to be proved in the next subsection. 

*.W = ^*.(^=.a). ..e. *.(. + i_f)=i-i*„(^(. + l-i),2-,). (51) 

*•<->- 2^*" (¥- ^) • '- *• h ^) ^ ¥*" ( 2^ (" ^) • ^ - ') ■ f'^^' 

These relations are an extension of the previously known functional conjugacy that relates s — A and s = —2. 

It is better notation to call the alternate series solutions ^* (x, s) and $* (x, s), instead of vj/^ (z) and $* (z). Then 
the "dual" parameter is 

s, = 2 - s , (53) 
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and the previous relations are more succinctly written as: 

**(a;,s) = — *o ( — ,s, ) , $, (x,s) = — $0 ( — ,s* ) • (54) 

That is to say, the series solution of Schroder's equation for map parameter 2 — s, about the trivial fixed point at 
a; = 0, is conjugate to the series solution of Schroder's equation for map parameter s, about the nontrivial fixed 
point at a;, = 1 — 1/s. Now, as previously noted, functional conjugation by g gives (PH) and (|l5|) . But then, 

g o so g^^ (u) = g [s — I + {2 — s) u) = *■ ^"Jy — h su , so 

gosog-^o-^g (w, s) = ^ _'^j + s^g {w, s) . (55) 

That is to say, the g-conjugated Schroder function obeys the inhomogeneous functional equation, 

^\~j'l + s^g (x, s) ^ M/g ((2 -s)xil-x),s) , (56) 

whereas ^* obeys a homogeneous equation, but with shifted position-dependent arguments. 

(2 - s) *, (x + ^— ^, s ) = ** (sx (1 - a;) + ^— ^, « ) • (57) 

To summarize the relations between the various functions: 

2 — s ( sx \ 

*o(a;,s) = ^— *J^— ^,2-sj , (58a) 

** (x, s) - ^ *o (^,2 - ^ , (58b) 

vl/g(x,s) = i^ + *, (^a; + |^,2-,s^ (58c) 

**(a;,s) = ^^ + 1'g (^a: + ^-^,2-sj (58d) 



1-s s /s-l + (2~s)a; 
23^+2-^*" 



*9(2^.5) = :^- + 7^-*o( ^ ^,s) , (58e) 



*o(x,s) = ^ + ^*J^^-^,s) , (58f) 



For the potential itself, functional conjugation leads to 



2-s\^, / s ( ( \ 



F(.,.)=^^— j r(_— (_.-(^l-jjJ.2-,J . (59) 

.,,,.(<i.4M^)V(^(._(._i)),_.)^ 

Perhaps it is worthwhile to work out the various functions explicitly for the cases which can be solved in closed- 
form, with ^0 and $0 = "^^^ as given in ([T7| . Consider the first and the last of these. The corresponding dual and 
conjugated functions are given by: 

**(a^,4) = -i*o(-2x,-2) = -^h7r-3arccos(^2.T + ij j , (61a) 



1 .\ ^/ . / 1 

a 



*4x,-2) = -2*o(-22;'4j =-2 arcsinW--a;l , (61b) 



*g(x,4) = ^ + vl/, (^x-|-2J =^-2farcsinjiQ-xjj , (61c) 

*3 (a;, -2) = ^ + ^, ( X - ^, 4 j === ^ - ^ (27r - 3 arccos (2a; - 1)) . (61d) 

We plot these, say for x such that the arguments of the inverse trigonometric functions are real, and the functions are 
principal- valued. 
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*o {x, 4) red, *, (x, 4) blue, and *g {x, 4) green. *o (x, -2) red, ** {x, -2) blue, and *g (x, -2) green. 



B. Alternate series solution theorem 

Consider the Poincare functional equation for the Schroder function inverse, rendered to facilitate expansion about 
the nontrivial fixed point "*" corresponding to x* = 1 — 1/s, 



$* ((2-s)z,s) = (2-s)$* (2,s)-s$^ (z,s) with $4z, s) = z + O (z^) , 



where 



Let 



X = I hz, Z = X -\ 1 

s s 



2 — s 

<&* {z,s) = (pizjs) 

s 



to change the functional equation to 

2-s 



0((2-s)z,s) 



(2-,sr 



</.(z,A)(l-0(z,A)) 



That is to say, 



sz 



(/)((2-s)z,s) = (2-s)(/)(z,s)(l-0(z,s)) with </)(z,s) = - 1" O (z^) 



2-s 



(62) 
(63) 

(64) 

(65) 
(66) 



Comparing this last relation to the original Poincare functional equation tailored to yield the expansion about the 
trivial fixed point "0" , 



^oisx,s) = s^oix,s){l- ^oix,s)) with $o (a;, s) = a; + O (a;^) , 



it follows that 



So then 



(t)iz,s) = $0 



2-s 



,2-s 



^ , X 2 — s^/sz^ \ ^ , , 2 — s^/sz^ 

$4z,s) = — — $0 ( ^3^,2-sj , $o(z,s) = — — $, ( ^— -,2-s 

and therefore the Schroder function, i.e. the inverse of the inverse function ^ = $^^, is 

T/ X 2 — s^/sz^ \ ^ , . 2 — s^/sz„ 

^4-,^) = ^*o(^^3^,2-sj , v,„(,,,) = ^v,,(^-_,2-s 

Hence $* (*, (z, s) ,s) — z ~ "ii^, ($, (z, s) , s). This proves the relationships between $,, *, and $o, ^o- 



(67) 



(68) 



(69) 



(70) 
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VII. APPENDIX B: CONSTRUCTING THE BRANCHES OF SCHRODER'S * FUNCTION 

We wish to imitate for ^ {x, s) what we did to find the branches of the potential in Section III of the text. We start 
with the functional equation, (|S]). Next we write y — sx (1 — x), with solutions x± = -^ il zt ^\ — 4j//s] . Now, we 
rename y — >■ a;, so the functional relation becomes 



*± {x, s)^s^\-[\± v/l^4a7s ) , s 



(71) 



One of these (^-) reproduces the original auxiliary function expanded about x — 0, only it does so more accurately 
for X — >■ s/4 as may be seen by numerical evaluation, while the other (vE'+) gives the auxiliary on another sheet of 
the function's Riemann surface. This is useful for determining the first switchback potential. The process may be 
repeated to get the other branches of the auxiliary as a sequence of functions, ^„. 

Here we are especially interested in s = 10/3, as we wish to confirm the order in which the switchback potentials 
are encountered by the evolving particle for this case. We find the following numerical results. 



Psi lO" 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 



Eight branches of * (a;, 10/3). 

It follows that the order of the potentials is as given in the text: Vb — ^ V^i — J> V2 — > T4^i — > W2 — > V3 — > Xi — > X2 — > • • • . 
As a practical matter, it is somewhat easier to compute the inverse function $ = ^^^ by combining series solution 
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methods with functional extensions. The functional equation for the inverse of Schroder's function is 

$ {sx, s) = s^ {x, s) (1 - $ {x, s)) . 
For s < 1 this quadratic equation may be solved for $ (x) in terms of $ (sx), 



*± (a;, s) = - 1 ± W 1 - -$ {sx, s) , 



(72) 



(73) 



and, upon iteration, all values of $ may be obtained from those values given accurately by the series about a; = 0. 
For s > 1, on the other hand, we first rescale x in ([7^ to write 



$(x,s)=.$(^,.)(l-$(J,s)) , 



(74) 



and from this equation, upon iteration, all values of $ may be obtained directly from those values given accurately 
by the series about a; = 0. 

As an example, consider again the case s — 10/3. We find: 



Phi 



0.8 



0.7 



0.6 



0.5 



0.4 



0.3 



0.2 



0.1 



0.0 




10 



20 



30 



40 



50 



50 



70 



80 



90 



100 
X 



Seventh order series approximation to $ (x, 10/3), in black, compared to results from using this series in the 1st 
(red), 2nd (blue), and 3rd (green) iterates of the functional extension ((74)) . 

Flipping this graph about the SW-NE diagonal, we obtain a graph for the multi- valued ^ = ^^^. From the third 
iterate of (1741) (as given in the graph by the green curve) the resulting curve for ^ (a;, 10/3) agrees with the previous 
graph obtained by solving Schroder's equation through use of the series about x — 0, combined with functional 
methods. 



VIII. APPENDIX C: THE s = 1 SERIES FOR THE POTENTIAL 



For this case the appropriate limit of (jlip gives 
V {x, 1) EE - lim ((In^ s) U {x, s)) 



= -x^ - 2x^ - 4x^ 



25 



215 



589 
l5" 



7813 
90 



xi"- 



60 481 
315 



x"- 



11821 
^8~' 



Oi 



„13^ 



(75) 
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More systematically, let 

/ oo N 

V {x, 1) = -xM 1 + ^ c„ x" 



Ci = 2 , C2 = 4 , C3 = 



25 



(76) 



and solve by iteration the functional equation that should be obeyed by V (x, 1), namely, 

V{x{l-x),l)^{l- 2xf V (x, 1) . 



(77) 



If the formal series ([751) is constructed to O (x^^) or so, it appears to have a radius of convergence of i? ~ 1/2. But 
remarkably, this is seen to be illusory to higher orders. Different behavior sets in around O (x"^"), where the successive 

|c„| " used in the limsup determination of R, p^. begin to grow linearly for n > 30, as shown here (purple curve, 
with wiggles). 




.l/r 



for s = 1 and 1 < ?i < 300. The blue horizontal is 2, while the sienna line is (2 "/^e '^"/^n!) 



l/r, 



That is to say, |c„| ~ L" exp (nlnn) for large n, with L : 
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5. 804 29 X 10 ■^. This immediately brings 



270 — V2eV2 

to mind Stirling's formula, n! ~ -\/27rn (^) , and in fact, a comparison of c„ with /„ = 2^"/^e~^"/^n! is striking, as 
shown in the Figure. (The most important feature in the Figure is agreement between the averaged slopes. A better 
overall fit to the c„, on average, is achieved by shifting the sienna line slightly to the left, for example by including 
an additional ^/n factor in /„.) 

All this is compelling evidence that ([75]) itself is not convergent, but rather an asymptotic series. Indeed, the /„ 
"fit" to c„ can be obtained from an asymptotic approximation of a simple integral: 



I{x)= I e 



1 



n=0 



%/2 /3x-2V2ei 

i^'^y = -— c'^p I — ^ — 



/2^ 



Ei( l,--\/2e^ 

X 



V Sr^ T,no-"/2„-3n/2 n 



x"2 



E 

n=0 



(78) 



Taking the Cauchy principal value for / (x) gives finite numerical results for all x. These results may be used to 
compute corrections to the polynomial approximation for V (x, 1) as constructed from truncating the series (j76p . We 
leave further analysis of this interesting but peculiar case to the reader. 
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